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Rates for detection of weakly-interacting massive-particle (WIMP) dark matter are usually carried 
out assuming the Milky Way halo is an isothermal sphere. However, it is possible that our halo is not 
precisely spherical; it may have some bulk rotation; and the radial profile may differ from that of an 
isothermal sphere. In this paper, we calculate detection rates in observationally consistent alternative 
halo models that produce the same halo contributions to the local and asymptotic rotation speeds 
to investigate the effects of theoretical uncertainty of the WIMP spatial and velocity distribution. 
We use self-consistent models to take into account the effects of various mass distributions on the 
local velocity distribution. The local halo density may be increased up to a factor of 2 by flattening 
or by an alternative radial profile (which may also decrease the density slightly). However, changes 
in the WIMP velocity distribution in these models produce only negligible changes in the WIMP 
detection rate. Reasonable bulk rotations lead to only an 0(10%) effect on event rates. We also 
show how the nuclear recoil spectrum in a direct-detection experiment could provide information 
on the shape and rotation of the halo. 
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I. INTRODUCTION 

Perhaps the most intriguing explanation for the dark 
matter in the Galactic halo is that it is composed of 
weakly-interacting massive particles (WIMPs) §• These 
particles typically have masses between 10 GeV and 
a few TeV and couple to ordinary matter only with 
electroweak-scale interactions. For example, the leading 
candidate WIMP is perhaps the neutralino, the lightest 
superpartner in supersymmetric extensions of the stan- 
dard model Several complementary efforts are cur- 
rently afoot to detect these halo dark-matter particles. 
For many WIMP candidates, the most promising av- 
enue is direct detection of the O(10 keV) recoil energy 
deposited in a low-background laboratory detector when 
a halo WIMP scatters from a nucleus therein B0]. An- 
other promising technique for many other WIMP can- 
didates is detection of the energetic neutrinos produced 
by annihilation of WIMPs which have been captured in 
the Sun and/or Earth ||. There are also efforts to de- 
tect anomalous cosmic-ray positrons, antiprotons, and 
gamma rays which may have been produced by WIMP 
annihilation in the Galactic halo (see Ref. for a review 
and further references). 

The predicted rates for all of these techniques depend 
on the mass and interactions of the WIMP. The rates for 
scattering from nuclei also depend on quantities such as 
quark densities in the nucleon and on nuclear form fac- 
tors. Considerable efforts have been made to survey the 
plausible parameter space for supersymmetric WIMPs. 
Furthermore, the sources of uncertainty in the predicted 
direct-detection and energetic-neutrino rates from, e.g., 
quark densities and nuclear form factors have been eval- 
uated and isolated. 

Of course, predictions depend on the spatial and ve- 
locity distribution of WIMPs in the halo. In most (all?) 
calculations of dark-matter detection rates, the halo is 
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assumed to be a cored isothermal sphere parameterized 
by a central (or alternatively, local) density and core ra- 
dius which are fit to the observationally inferred halo 
contribution to the rotation curve. In this model, the 
velocity distribution is Maxwcll-Boltzmann with a veloc- 
ity dispersion determined by the rotation speed at large 
radii. Observational uncertainties in the rotation curve 
and in the disk and bulge contributions lead roughly to 
a factor-of-two uncertainty in the local dark-matter den- 
sity. Assuming an isothermal sphere, one finds the local 
dark-matter density po = 0.2 — 0.4 GeV cm -3 and a ve- 
locity dispersion v = 270 ± 70 km sec^ 1 . 

In addition to these uncertainties from the rotation 
curve and disk mass distribution, deviations from the 
standard nonrotating isothermal spherical halo are also 
plausible, if not probable. Essentially all the empirical 
information we have on the halo is provided by the ro- 
tation curve. To a first approximation, almost any halo 
mass distribution which gives rise to a flat rotation curve 
is acceptable. Although there are some arguments that 
the halo must be more diffuse than the disk ||, there is 
no reason why it should be perfectly spherical. In fact, 
there is ample evidence that the halos of several external 
spiral galaxies are flattened by roughly a factor of two 
and now some evidence that the Milky Way halo is 
similarly flattened ||. The dominant effect of flattening 
on the detection rate is through the local dark-matter 
density j^]. However, flattening may also affect detec- 
tion rates through the velocity distribution, which has 
not been taken into account. 

Bulk rotation can also affect the velocity distribution 
of WIMPs seen at the Earth. Again, the rotation curve 
is determined by the halo mass distribution and is insen- 
sitive to its velocity distribution. Therefore, there is no 
empirical evidence to rule out a halo with some bulk rota- 
tion. Although there are theoretical arguments against a 
rotation-dominated velocity distribution, there are also 



1 



reasons to expect the halo to have some bulk rotation 

There may also be theoretical uncertainties in the halo 
radial profile. The functional form for the radial profile 
commonly assumed is in fact a phenomenological model 
which produces a linearly rising rotation curve at small 
radii and a flat rotation curve at large radii. There are 
other radial profiles which will satisfy these requirements 
and produce the same rotation speeds at the Galactocen- 
tric and large radii to which the models are fit. 

In this paper, we investigate uncertainties in the 
WIMP detection rate which arise from imprecise knowl- 
edge of the spatial and velocity distribution of dark- 
matter particles. To do so, we use a class of self- 
consistent models for a flattened and /or rotating halo 
which have been developed by Evans 113] , and consider 
several plausible spherical distributions. All the models 
we consider produce the same halo contribution to the lo- 
cal and asymptotic rotation speeds. Some of our models 
may appear to be extreme (in terms of flattening, bulk 
rotation, or central density) to some Galactic-dynamics 
experts; however, our primary aim is to provide a conser- 
vative estimate of the uncertainty in dark-matter detec- 
tion rates from uncertainties in the halo distribution, and 
the models we consider span a range of observationally 
plausible — though not necessarily theoretically favored — 
models. 

We find that flattening and/or changes to the radial 
profile may increase the density by roughly a factor 
of two. However, either departure from the canonical 
isothermal sphere has a negligible effect on the velocity- 
distribution dependence of the event rate. The bulk ro- 
tations which may arise in realistic galaxy-formation sce- 
narios will have no more than a 10% effect on detection 
rates. 

In the next Section, we review the procedure for cal- 
culating detection rates. In Section III, we review the 
distribution functions for the Evans models which we use 
to investigate the effects of flattening and bulk rotation. 
Results for the effects of flattening on the local WIMP 
velocity distribution, density, and total and differential 
detection rates are provided in Section IV. We also pro- 
pose here that the measured differential recoil-energy dis- 
tribution (in case of detection) could be used to constrain 
the bulk rotation and flattening of the halo. In Section V 
we investigate the effects of uncertainties in the halo ra- 
dial profile in spherical models on dark-matter detection 
rates. In Section VI we summarize and make some con- 
cluding remarks. We also discuss how rates for indirect 
detection of WIMPs will be affected in these alternative 
halo models. 



II. CALCULATING DIRECT-DETECTION 
RATES 

One can write the differential rate for direct WIMP 
detection Jl| as 



dQ 2m x mt J„ mm v 

where do is the cross-section (at zero momentum trans- 
fer); po is the local dark matter density; m r is the reduced 
mass mNm x (mN+Tn x )~ 1 , where tun is the mass of a tar- 
get nucleus and m x is the WIMP mass; Q — |q| 2 /2mjv is 
the deposited energy, where q is the momentum transfer; 
F(Q) is a nuclear form factor; fi(v) is the distribution of 
WIMP speeds relative to the detector (normalized to 1); 
and Vmi n = [(QTOAr)/(2TO 2 )] 1 / 2 . Defining the dimension- 
less quantity, 

T(Q) = ^vo r ^dv, (2.2) 

and taking F(Q) = exp(— Q/2Qq), the differential detec- 
tion rate can be written as 

dR = f Po «o _ \ eM _ Q/Qo)T{Q) . (23) 
dQ \y/7TVQm x m^ J 

i.e., the density times a velocity-dependent term. The 
total event rate can be determined by integrating over 
all detectable energies: 

r°° dR 

where £t is the threshold energy for the detector. 

III. HALO MODELS 

To study the effects of flattening and bulk rotation on 
detection rates, we use Evans's family of analytic axisym- 
metric distribution functions (DFs) |T^ ], 

F(E, I 2 ) = [AL 2 Z + B]exp(AE/v 2 ) + Cexp(2£/« 2 ), 

(3.1) 

with 

/2\ 5 / 2 (l_a p _{2 \ 1/2 R 2 V-l 

(3.2) 

where E is the binding energy, L z is the azimuthal com- 
ponent of angular momentum, vq is the circular speed at 
large radii, R c is the core radius, and q is the flattening 
parameter, ranging from q = 1 for a cored, spherical halo 
to q = l/y/2 ss 0.707 for the most flattened non-negative 
DF |12| ). These models elegantly reproduce Binney's po- 
tential and corresponding density |13| , 

^(R,z) = -±v 2 \ogfRl + R 2 + ^\ (3.3) 
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where R is the radial distance and z is the vertical dis- 
tance above the disk. These are suitable for describing 
the halo since they produce rotation curves which rise 
linearly at small radii and are flat at large radii. 

In this calculation, we take the dark-matter contribu- 
tion to the local circular velocity to be v c (Rq) = 170 
km sec -1 (which we get from a local rotation speed of 
220 km sec -1 and a disk contribution of 140 km sec -1 ), 
vq = Vqc = 220 km sec -1 , a Galactocentric radius 
Ro = 8.5 kpc, and z = kpc. The core radius R c is 
obtained by noting that in the plane z = 0, 



# = v 2 R 2 
dR R 2 + R 2 ' 



Therefore, for all q, the core radius is 



R c — Ro 



Vc(Ro) 



1/2 



7 kpc. 



(3.5) 



(3.6) 



[We have checked that our conclusions on the effects of 
flattening are unchanged if wc adopt other plausible val- 
ues for Vqo, v c (Ro), and Ro.] 

The isopotential contours for these models are ellip- 
soidal with (short-to- long) axis ratios q [c.f., Eq. (3. 3)]. 
Fig. |l| shows isodensity contours for q = 1, 0.85, and 
l/y/2 « 0.707 (for R c = 7 kpc). The isodensity con- 
tours are not ellipsoidal. For small radii, they are close 
to spherical, and they become more flattened for larger 
radii. Inspection of Fig. [I] shows that (for R c = 7 kpc) 
the short-long axis ratio for the isodensity contours is 
roughly 1:2 for q ~ 0.707 for radii comparable to our 
Galactocentric radius. 

The DFs above have no bulk rotation. However, a 
family of DFs with bulk rotations can be constructed by 
considering linear combinations, 

G(E,L Z ) = aF+(E,L z ) + (l-a)F_(E,L z ), (3.7) 

of DFs 



+ v z ' \0, < 0, 



F-{E,L 2 Z ) = 



0, > 0; 

F(E,L% « <O, 



(3.8) 



(3.9) 



with only positive or negative azimuthal-velocity compo- 
nents Vcf,. These models have the same spatial distribu- 
tions as the nonrotating models F(E,L Z ). The parame- 
ter a ranges from 1 (for maximal corotation) to 0.5 (the 
model with no net rotation) to (maximal counterrota- 
tion), and is related to the dimcnsionlcss spin parameter 
A usually used to quantify galactic angular momenta by 
A = 0.36|a- 0.51. 
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FIG. 1. Halo isodensity contours for the Evans models 
q = 1, 0.85, 0.707, where the z and R axes are in kpc.) 
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The DFs discussed so far specify the velocity distribu- 
tion in the Galactic rest frame. However, the solar system 
moves with respect to this frame with a velocity v s = 220 
km sec -1 . Therefore, the DF F s (vr, v z , v^) with respect 
to the Earth can be obtained from the rest-frame DF F 
by substituting — * + v s . 



A. No Net Rotation 

For these models, the distribution function is even 
in the variable w^; there are as many particles circling 
around clockwise as there are counterclockwise. 

Substituting the binding energy, 

E = -\v 2 - \vl log (r\ + R 2 + z 2 /q 2 ^j , (3.10) 

into the distribution function F(E, L 2 ) and transforming 
to the Sun's rest frame yields 

F S (E : L 2 Z ) = [AR 2 { V(j> +v s f + B] 

e^{-^(v 2 +v 2 + (v^+v s ) 2 )} 



C 



{R? c + R 2 + z 2 /q 2 ) 2 
e^{~^(v 2 r +v 2 + ( V4l + v s ) 2 )} 
{R 2 C +R 2 + z 2 /q 2 ) ' 



(3.11) 



Since v 2 = v 2 R -\-v 2 + v 2 1 , one can simplify this to depend 
only on v and the angle a between the velocity and the 
azimuthal direction. Plugging in for the local coordinates 
(i?, z) = (Rq,0), one obtains the more convenient form, 



f(v,a) = [ARl(v cos a + v s ) 2 + B] 

exp{— -^i{v 2 + 2w s ucosa + v 2 )} 

x ° (Rl + Rl) 2 

exp{— -jz(v 2 + 2u s t!Cosa + v 2 )} 

{Rl + Rl) ' 



(3.12) 
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where the q dependence is still implicit in the coefficients. 
Therefore, the local speed distribution function needed 
for calculation of the dark-matter detection rate is 



fi(v) 



fn f( v : a ) U 2 sm a dd 

Io° Jo f( v ' a ) 1,2 sin a ^ a dv 



(3.13) 



The top panel in Fig. g shows the speed distributions 
fi(v) for the nonrotating halo for q — 1, 0.85, and 0.707. 
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B. Maximally Corotating and Counterrotating 

The calculation of the speed distribution for a rotating 
halo proceeds in the same fashion. However, for the max- 
imally corotating case, the DF in the Sun's rest frame is 

F„<E,ri) = »*;;»;; (3.i4) 



FIG. 2. Local speed distributions fi(v) for nonrotating, 
maximally corotating, and maximally counterrotating mod- 
els with q = 1 (solid curves), q = 0.85 (dashed curves), and 
q = 0.707 (dotted curves) and with v s = 220 km/s. 
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and for the maximally counterrotating model, the DF in 
the Sun's rest frame is 



0, 

Fs(E,L 2 z ), 



(3.15) 



The middle and bottom panels in Fig. g show the 
speed distributions fi(v) for the maximally-corotating 
and counterrotating models, respectively, again for q = 1, 
0.85, and 0.707. Note that there are no particles with 
v < v s for the maximally counterrotating model. Also, 
the steep rise in fi(v) near v = for the q = 0.707 
corotating model arises because there are more particles 
in nearly circular orbits with velocities v s — nearer to 
in our frame — in this model than in the q = 1 model. 
The maximally rotating models have a spin parameter 
A = 0.18 which is significantly larger than the spin pa- 
rameters A ~ 0.05 expected from galaxy-formation mod- 
els |llj. Therefore, realistic speed distributions should lie 
somewhere between these two and closer to that for the 
nonrotating model. 



IV. TOTAL AND DIFFERENTIAL 
DIRECT-DETECTION RATES 

Fig. U shows the differential detection rates dR/dQ 
for spherical and flattened nonrotating and maximally 
corotating and counterrotating models. It is seen that 
flattening has a weak effect on the predicted differential- 
detection rate. Bulk rotation (especially counterrotation) 
has a somewhat stronger effect on the differential rates. 
Therefore, the shape of the nuclear recoil spectrum could 
provide information on whether the halo is rotating or 
not, and this could be useful for constraining galaxy- 
formation models. 

Fig. H shows the total detection rate (assuming no 
thresholds) for nonrotating and maximally corotating 
and counterrotating models as a function of the flattening 
parameter q. The detection rate increases roughly as q~ x 
independent of the rotation. The larger incident WIMP 
velocities in counterrotating models leads to a stronger 
form-factor suppression. This is the leading factor in 
accounting for the decrease in the event rate in coun- 
terrotating models and vice versa for corotating models. 
Maximal rotation can change the event rates by roughly 
30%. However, the spin parameters expected on theoret- 
ical grounds are generally smaller than a third of that for 
our maximally rotating halos. Therefore, the most plau- 
sible values for the bulk rotation should yield detection 
rates within 10% of those for the canonical nonrotating 
model. 

The q depend enc e of the local halo density can be de- 
rived from Eq. ( |3.4| ), and Fig. || shows that it scales very 
nearly as q~ x . Fig. || shows the detection rate scaled 
by the local halo density as a function of q. These two 
Figures illustrate that the change in the velocity distri- 
bution from flattening has essentially no effect on the 
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FIG. 3. Differential detection rates for Evans's models 
q — 1 (top) and q = 0.707 (bottom) with no net rotation 
(solid), maximal corotation (dashed), and maximal counter- 
rotation (dotted). 
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FIG. 4. Total detection rate as a function of halo flattening 
q for nonrotating (solid), maximally corotating (dashed), and 
maximally counterrotating (dotted) halos. 
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dark-matter detection rate. Heuristically, halo particles 
move in the same gravitational potential as the Sun, and 
the velocity dispersion of any species is fixed by the po- 
tential. Our calculations contradict the claims of Cowsik 
et al. IT3] and verify the arguments of Refs. Jig] . 

For these calculations, we have used a WIMP with 




only scalar interactions of mass m x — 100 GeV and 
(To = 4 x 10~ 36 cm 2 and a Germanium target nucleus. 
We have checked that our conclusions do not change if 
we use a different WIMP mass and/or target nucleus. 
We have also checked that this conclusion is independent 
of the details of the assumed rotation curve: The veloc- 
ity dispersion is essentially independent of the flattening 
in models where the halo contribution to the local rota- 
tion curve is higher or lower than that which we have 
used here, either because of different measured rotation 
speeds, or because of a different disk/bulge contribution. 
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FIG. 5. Local halo density as a function of the flattening 
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Let us now consider the effect of possible variation 
in the radial profile in spherical halo models. Heuris- 
tically, galaxy formation results in a cored isothermal 
halo through the process of violent relaxation. How- 
ever, there will realistically be some deviations from this 
simple physical picture for halo formation. For exam- 
ple, the collapse of baryonic matter in the Galaxy might 
affect this process. Empirically, the evidence for flat- 
tened spiral-galaxy halos suggests some departure from 
the simple picture. Therefore, even if we consider only 
spherical halo distributions, there is still some latitude in 
our choice of the precise form for the radial profile of the 
halo. 

An empirically plausible radial profile for a spherical 
Galactic halo is constrained by its contribution to the 
Galactic rotation curve. Therefore, it should approach a 
constant near the core so it gives rise to a linearly rising 
rotation curve at small radii, and it should fall as r~ 2 at 
large radii to provide a flat rotation curve. The canonical 
profile usually used for dark-matter calculations is the 
so-called "isothermal" sphere (actually, the radial profile 
of the true cored isothermal sphere cannot be written 
analytically; see Rcf. p. 229), 



p(r) 



(5.1) 



FIG. 6. Velocity dependence (i.e., the detection rate scaled 
by the local halo density) of the total detection rate as a 
function of the flattening q for nonrotating (solid), maximally 
corotating (dashed), and maximally counterrotating (dotted) 
halos. 



where r c is a core radius which is fit to the halo con- 
tribution to the local rotation speed (and r is now the 
spherical radial coordinate, r 2 = R 2 + z 2 ). Of course, 
the radial profile of the spherical Evans model, 



p(r) 



3r 2 



AttG (r 2 +r 2 ) 



2\2 ' 



(5.2) 



also has the desired properties. Yet another analytic form 
which might be empirically acceptable is 
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FIG. 7. The radial profile of the three spherical halo mod- 
els discussed in the text. The solid, short-dash, and long-dash 
curves are for the canonical isothermal [Eq. (5.1)], spherical 
Evans [Eq. (5.2)], and alternative isothermal [Eq. (5.3)] mod- 
els, respectively. The dotted curve is for the non-increasing 
radial profile that gives the largest local density. 
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Keep in mind that the core radius r c for each model must 
be fit to the rotation curve, and r c for each model will 
be different. Suppose, as we did before, that the local 
rotation speed is 220 km sec -1 and the disk contribution 
is 137 km sec -1 . Then the local halo contribution to the 
rotation curve is 170 km sec -1 which leads to core radii 
r c = 7 kpc (as before) for the Evans model, r c = 2.8 kpc 
for the canonical isother mal sphere, and r c — 0.9 kpc for 
the alternative in Eq. (5.3). The rotation curves and 
radial profiles for these three models are plotted in Figs, 
[j] and U respectively. Th e s olid , short- dash , and long-dash 
curves are for Eqs. (5.1), (5.2), and (5.3), respectively. 

In principle, the radial profile can be fixed by the halo 
contribution to the rotation curve. However, measure- 
ment of the Galactic rotation curve is notoriously diffi- 
cult, especially near the interior. Furthermore, the disk 
contribution to the rotation curve must be known to in- 
fer the halo contribution, and precise determination of 
the disk contribution is also difficult. Therefore, there 
will be significant uncertainties in any reconstruction of 
the halo contribution to the rotation curve from observa- 
tional data. 

Using our canonical values for the halo contributions to 
the rotation speed, we find local halo densities of 0.43,[] 
0.51, and 0.38 GeV cm -3 for the isothermal, Evans, and 
alternative models respectively. Therefore, although the 
central density of these three models may differ consid- 
erably (c.f. Fig. 0), the requirements that each yield the 
same halo contribution to the local rotation speed and 
the same asymptotic rotation speed constrain the local 
halo density in these models to 20%. 

One could contemplate a profile with a smaller local 
density with a higher central density. However, a profile 
with a central density much greater than that in our al- 
ternative model will have a core density comparable to 
the Bulge density (approximately 50 GeV cm -3 and 
will therefore contradict observed Bulge dynamics. It is 
therefore unlikely that the local halo density can be re- 
duced while maintaining the same halo contribution to 
the local and asymptotic rotation speeds. Contrariwise, 
one could consider a model with a larger local density and 
smaller core radius which still gives the same contribution 
to the local halo speed. Any physically reasonable radial 
profile should be monotonically decreasing with radius. 
The limiting case (a density which is constant interior 
to our Galactocentric radius; the dotted curve in Figs. |7] 
and ||) yields a density 1.4 GeV cm -3 [v c (rQ)/v oa ] 2 , which 
results in 0.8 GeV cm -3 for a local halo rotation speed 
of 170 km sec" 1 . Therefore, a local halo density roughly 



FIG. 8. Rotation curves for the spherical models shown in 
Fig. 7. 



*This corrects the value of 0.35 GeV cm 3 given in Section 
2.4 of Ref. §. 
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twice that obtained from the canonical model is conceiv- 
able (although perhaps somewhat artificial as indicated 
in Fig. 0), and a local halo density 0(10%) smaller than 
the canonical value is also possible. 

We have evaluated numerically the direct-detection 
rate using the DF, Eq. (fo.lf ), for the spherical Evans 
model and the Maxwell-Boltzmann velocity distribution 
for the isothermal sphere. We find that the detection rate 
with the spherical Evans model is roughly 15% larger 
than that in the isothermal model. Therefore, the dif- 
ference in detection rates can be attributed primarily to 
the difference in the local halo density and only secondar- 
ily to the differences in the velocity distribution. Once 
again — as in the case of flattening — we find that different 
radial profiles lead to roughly the same velocity disper- 
sions as long as both profiles are fit to the same halo 
rotation speed. 

VI. CONCLUSION 

Predictions for WIMP detection rates are almost al- 
ways carried out assuming the dark-matter distribution 
to be an isothermal sphere. When fit to reasonable values 
of the halo contribution to the local and asymptotic ro- 
tation speeds, the canonical isothermal halo gives a local 
halo density 0.25-0.5 GeV cm~ 3 . Its velocity distribu- 
tion is Maxwell-Boltzmann with a velocity distribution 
fixed by the asymptotic rotation speed. 

However, virtually all the empirical constraints to the 
halo come from its observationally inferred contribution 
to the Galactic rotation curve. These (still rather poorly 
determined) data are supplemented by some qualitative 
theoretical notions about the halo: i.e., that it should be 
more diffuse than the disk and monotonically decreas- 
ing with Galactocentric radius. Many halo distributions 
can satisfy these observational and theoretical constraints 
and still produce the same local and asymptotic rotation 
speeds. 

In this paper, we have calculated WIMP direct- 
detection rates in several plausible alternatives to the 
canonical model. We find that if the halo is flattened 
with an isopotential axial ratio q, the direct-detection 
rate will increase by roughly q^ 1 . This increase is due 
primarily to the effect of flattening on the local halo den- 
sity, which also increases as q^ 1 . We have used a self- 
consistent distribution function for a flattened halo to 
verify that the effects of flattening on the velocity dis- 
tribution have virtually no effect on the detection rate. 
Local stellar kinematics and the thickness of gas layers 
suggest that halo isodensity contours may be flattened 
by up to a factor of 2 || corresponding to q = 0.707 — 1 
for the Evans models, which would suggest that flatten- 
ing might increase the local halo density, and therefore 
direct-detection rates, by a factor of 1.4. However, the 
heuristic argument that flattening should affect the de- 
tection rate primarily through its effect on the density, 



and only secondarily through its effect on the velocity 
distribution should also apply to a halo with ellipsoidal 
isodensity (rather than isopotential) contours. In such 
models, the local density is increased by a factor near 2 
for a flattening near 2 || . 

There are no empirical constraints to the bulk rotation 
of the halo. A maximally corotating or counterrotat- 
ing halo could increase or decrease the detection rate by 
40%. However, galaxy-formation scenarios generally pre- 
dict bulk rotations no more than 0.3 of maximal. There- 
fore, we do not expect bulk rotation to change the pre- 
dicted event rates by more than 10%. Although simple 
galaxy-formation models suggest that a halo would coro- 
tate if it rotated at all, the existence of counterrotating 
disks JTtJ suggests that a counterrotating halo might also 
be plausible. 

We found that in spherical models, the local density 
could be increased by up to a factor of two and decreased 
slightly with different radial profiles that still give the 
same local and asymptotic halo rotation speeds. In this 
work, we focussed on halos with axial symmetry, but it is 
possible that the halo may deviate somewhat from axial 
symmetry. However, detection rates in reasonable triax- 
ial models also generally fall within a factor of two of the 
canonical detection rates Q . 

We restricted our analysis to direct detection. How- 
ever, similar conclusions should apply to rates for indirect 
detection of WIMPs via observation of energetic neutri- 
nos from WIMP annihilation in the Sun and/or Earth. 
Like direct-detection rates, these rates are controlled pri- 
marily by the local halo density. Since the velocity dis- 
persion is fixed to a large extent by the local and asymp- 
totic rotation speeds, indirect-detection rates should not 
be affected by their dependence on the velocity distribu- 
tion. 

On the other hand, plausible deviations from the 
canonical isothermal sphere can lead to dramatically 
different fluxes of anomalous cosmic-ray antiprotons, 
positrons, and gamma rays from WIMP annihilation in 
the halo. These fluxes are determined by an integral of 
the square of the density over the entire halo. Although 
the local halo density does not differ too much in alterna- 
tive models, the core density can differ dramatically. In 
particular, Fig. [?] shows that the central density can be 
increased perhaps by an order of magnitude over that in 
the canonical model. If so, then the flux of gamma rays 
from WIMP annihilation in the Galactic center would be 
increased by a factor of 100 over the fluxes predicted in 
canonical models. 

There is also the possibility that if WIMPs are de- 
tected, the nuclear recoil spectrum might tell us about 
the structure of the halo. Fig. || shows how the recoil 
spectrum could be used to constrain the rotation of the 
halo. We have also investigated the magnitude of annual 
modulations in the event rate due to the Earth's orbital 
motion around the Sun. We found a maximally corotat- 
ing halo could increase the annual modulation by a fac- 
tor of 2, implying an increase in modulation amplitude 
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of 0(30%) for models with more realistic corotation. 
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